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Abstract: 

We present if-band observations of the elliptical galaxy NGC 4696, the brightest member of the Cen- 
taurus cluster of galaxies. We have measured its light profile, using a two-dimensional fitting algorithm, 
out to a radius of 180 arcsec (37 h^^ kpc) . The profile is well described by a de Vaucouleurs law, with 
an effective radius of 35.3 ±1.0 kpc. There is no need for the extra free parameter allowed by a 
Sersic law. Allowing for a variation of 0.3% in the sky level, the profile obtained is compatible with 
data from 2MASS. The profile shows no sign of either a truncation or an extended halo. 
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1 Introduction 

Surface brightness profiles of elliptical galaxies are 
norm ally well described by either a de Vauc ouleurs 
law (|de Vaucouleurs I Il948|) or a Sersic law ijSersic I 
1968). However, the situation can be more com- 
plicated when we consider a galaxy embedded in a 
cluster. In this case, the outer parts of the profile 
may be modified due to tidal interactions with neigh- 
bouring galaxies, giving a 'truncated' profile. A dif- 
ferent case is often encountered in elliptical galax- 
ies that are the brightest cluster member (BCM) 
and are located at the cluster cent re, many of which 
have an extended envelope or halo ijg chombert Ill987l : 
iGonzalez. Zabludoff. fc Zaritskv II2005T) . These galax- 
ies are known as 'central dominant' (cD) galaxies. The 
extended envelope is observed as a deviation of the 
profile from a de Vaucouleurs or Sersic law at large 
radii, and appears to correspond to light from stars 
orbiting in the cluster potential and not bound to th e 
central galaxy (|Bahcall Ill977l : I Gonzalez et al.l feOQS'l. 
Simulations indicate that the stars in the envelope 
may, on average, be ol der than the stars in the galaxy 
ijMurante et al. Il2004h . 

NGC 4696 is the BCM of the Centaurus cluster, ft 
is an elliptical galaxy located near the centre of the 
cluster, ft is sometimes described as a cD galaxy, al- 
though no evidence has been found up to now of an 
extended halo. 

Previous measurements of the light profile of 
NGC 4696 at op tical wavelengths include th ose of 
ISchombert 1 J1987D and lJerien fe Dressier I (Il997fl. both 
obtain ed from photographic plates. I.TerienfcTjressler I 
l|l997l 'l measured the light profile out to a radius of 
~8 arcsec (~ 1 6fo7n kpc), while the profile obtained 
bv ISchombert I (|1987i ) reaches much further, out to a 
radius of ~9 arcmin (~ llOTlfg kpc). Neither observa- 
tion showed signs of an extended envelope. The profile 
from Sch ombert I (E?87) followed an r* law (de Vau- 
couleurs profile) out to a radius of «5 arcmin, and a 



truncation or 'distention' of the profile was found for 
larger radii. 

Measurements of the NIR profile were ob- 
tained by the Two Micron All Sky Survey 
(2MASS source 2MASX J12484927-4118399 
in the All-Sky Extended Source Catalog; 
http://www.ipac.caltech.edu/2mass/). Light profiles 
were obtained in three bands (J, H and K) for 
radii up to ~3 arcmin. This corresponds to surface 
brightnesses fin < 24 mag arcsec -2 in the H band. 

The Intra-Cluster Medium (ICM) of the Centau- 
rus cluster has been studied using X -ray emission by 
the R OSAT and ASCA satellites. lAllen fc Fabian I 
l|l994l '). using data from ROSAT, found a correla- 
tion between the structure — ellipticity and centroid — 
of the X-ray emission and optical observations 
flSparks. MacchettQ 1 _fc L _G_olorr^ekJl989T) of the central 
galaxy NGC 4696. Ilkebe et al. I lll999ft used data from 
both satellites, and described the X-ray emitting gas 
in the central region (r < 3 arcmin) according to a two- 
phase model. They identified the cool phase with the 
inter-stellar medium (ISM) of NGC 4696, and the hot 
phase with the ICM of the Centaurus cluster. These 
relationships between the central galaxy and the ICM 
could be an indication of the presence of an extended 
halo. 

The previous optical and NIR work pose a 
quandary in that the form of the profiles measured 
range from de Vaucouleurs to exponential. Thus the 
aim of the present work was to observe the light pro- 
file of NGC 4696 in the near-infrared (NIR) out to large 
radii, study the functional form of the profile (de Vau- 
couleurs or Sersic law) , measure the parameters of the 
profile, and seek evidence for either a truncation due 
to tidal interactions or an extended envelope following 
the cluster potential. The main advantage of the NIR 
over optical studies is that the light is dominated by 
old stars, and thus the NIR profile will better trace the 
bulk of the stellar mass of the galaxy. Furthermore the 
extinction due to dust is smaller in the NIR than at 



2 



Publications of the Astronomical Society of Australia 



optical wavelengths. 

Another motivation for this work is the ap- 
plication to studying radial colour distributions, 
through coupling with optical data, to understand 
the distribution of stellar populations throughout 
BC Ms. Through comparisons with simul ations 
fe.g. ISommer-Larsen. Romeo fc Portinari II2005T) . such 
studies would help to constrain BCM and cluster for- 
mation theori es and the enrichme nt of the intra-cluster 
medium fe.g. lLin fc Mohr 1120041) . 

We describe our observational data in Section [5] 
and explain the reductions and calibrations in Sec- 
tion El and the method used to fit the profile models in 
Section^] Section|K]presents the results and Section|S] 
summarises our conclusions. 

Throughout the paper we use a Hubble constant of 
Ho = 70 /i7o kms -1 Mpc -1 . We have used a redshift 
of z = 0.00987 ± 0.00005 llde Vaucouleurs et al. II1991I . 
which assumes no peculiar velocity) which equates to a 
distance of 42.3±0.2 h^ l Mpc; the angular scale at this 
distance is 0.205 ± 0.001 kpc arcsec - . Distances de- 
rived from surface brightness fluctuations, taking into 
account the peculiar motions of the Cen30 and Cen45 
components of the Centaurus cluster, and the Hydra 
cluster, give a distance of 42.5 ± 3.2 Mpc, consistent 
with the value used here (iMieske. Hilker. fc Infantel 
I2005TI . 

2 Observations 

Observations were made on 2005 February 16, using 
the AAO Infrared Imager & Spectrograph (IRIS2) at 
the 3.9m Anglo- Australian Telescope (AAT). Condi- 
tions during the night were photometric, but the see- 
ing was 2.5 arcsec. Data were obtained using the H 
filter. The detector in IRIS2 is a 1024x1024 Rockwell 
HAWAII- 1 HgCdTe infrared detector. The detector 
scale is 0.4486 a rcsec pixel" 1 , giving a field of view of 
7.7x7.7 arcmin jTinnev et al. 1120041) . 

We obtained a total of 99 frames of the field around 
NGC 4696 using a 9-point dither. Each of these frames 
is an average of 6 exposures of 10s each. To measure 
the light profile to the full extent allowed by the de- 
tector, NGC 4696 was placed in the south-east corner 
of the field of view. We also obtained observations 
of three photometric standard sta rs (UKIRT Faint 
Stand ards FS129, FS132 and FS135: lHawarden et al. 1 
2001), taking five 10s exposures in each case. In order 
to model the point spread function (PSF) , a bright star 
at a similar airmass to the science field was observed 
over 93 2.5s exposures. As explained in Sectional ex- 
tra data was needed for the fiat-fielding process. This 
was obtained from other //-band observations taken 
on February 11 during the same run. 

3 Reduction and Calibration 

The data reduction was carried out using the 
IRAF jhttp://iraf.noao.edu/| l and ORAC-DR 
(http://www.oracdr.org/ 1 software. The steps fol- 
lowed were: subtraction of dark current, flat-fielding, 



registering to obtain the correct astrometric offsets, 
and combination of the 99 frames. 

In the first step, the dark current in each pixel was 
subtracted from each of the raw frames. This value 
was obtained from a dark frame taken using the same 
exposure time (10s) as the science frames. The typical 
value of the dark current was ~3 ADU (~15 electrons) . 

To take into account the difference in response be- 
tween pixels in the detector, it was necessary to flat- 
field the frames. It was not possible to obtain a flat 
field directly from our science field, as the large extent 
of NGC 4696 meant that any flat field obtained from 
the science field images would retain some contribution 
from the galaxy. 

The flat field used was therefore derived from H- 
band observations made on February 11, during the 
same IRIS2 observing run. These were used instead 
of other observations from February 16 for a variety 
of reasons. In the case of the standard photometric 
stars, other very bright stars in the field caused arti- 
facts in the flat field. In the case of the 'PSF star', the 
accuracy of the resulting flat field was limited by the 
short exposure time. Using all the 10s exposed frames 
from February 11 that did not contain bright stars, we 
could construct a flat field with a total exposure time 
of about 58 min. 

The flat field from these frames was obtained us- 
ing ORAC-DR. For each of the fields, the frames were 
median-combined to obtain a first estimate for the flat. 
Then the original frames were divided by this estimate, 
objects were detected and masked, and the frames were 
combined again to obtain the final flat for the field. 
The flat used was an average of these flats. The dark- 
subtracted frames were divided by this flat before be- 
ing registered and combined. 

Before combining the 99 frames, registration was 
carried out to find the true offsets between the frames. 
This step was done using the xregister task in IRAF, 
which performs a cross-correlation between objects 
in the field. Knowing these offsets, the IRAF task 
combine was used to sigma-clip and median-combine 
the frames, obtaining the final image. The base sky- 
level is not subtracted, but rather an additive zero 
point is applied to each frame to account for differ- 
ences in sky level. This zero point was calculated as 
the mode of the same physical region of sky in all the 
frames. A 'sigma image' containing the standard devi- 
ation of the values corresponding to each pixel in the 
final image was also derived. 

The final image obtained using this reduction pro- 
cedure is shown in Figure Q Large-scale (of the order 
of arcmin) variations in the resulting background that 
do not correspond to the galaxy's light were observed 
at the level of 0.4% of the mean sky background, cor- 
responding to an uncertainty of ~ 0.004 mag arcsec -2 . 
These variations impose the effective limit on the ra- 
dius to which we were able to accurately measure the 
surface brightness profile. 

The photometric zero point for our observations 
was obtained from the photometric standards we ob- 
served. Our instrumental magnitudes for these stars 
were compared to their H magnitude in the 2MASS 
catalog. The zero point obtained has an associated 
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Figure 1: Final image used for the analysis. North 
is to the top and east is to the left. The field of 
view is 7.7 x 7.7arcmin. NGC4696 is located near 
the south-east (bottom left) corner. 



error of ±0.03 mag. 

4 Data Analysis 

Two further elements were needed prior to fitting the 
galaxy's profile: an estimate of the PSF in the observa- 
tions, and a mask for the other objects in the field. The 
profile could then be studie d using the two-dim ensional 
fitting algorithm GALFIT ijPeng et al. 112002ft . 

The PSF was obtained from the profile of a bright 
star observed immediately after NGC 4696 at a similar 
airmass to the science field, and hence with compara- 
ble seeing. A region of the image centered on the PSF 
star and covering f50xf50 pixels (~ 67 x 67 arcsec) 
was used directly to model the PSF. This PSF was 
convolved with the different galaxy profile models be- 
fore comparing them with the data during the fitting 
process. 

In order to test the robustness of the fitting pro- 
cess to the exact PSF, the fits were repeated with ar- 
tificial PSF images with varying FWHM, and Moffat 
and Gaussian profiles. Big differences were found be- 
tween a Moffat and Gaussian profile, but the FWHM 
had only a small effect (« 2.5%) on the fitted value 
of r e . Since we used an image of a star for the PSF, 
the shape of the profile will be accurately determined, 
and thus our fits should be robust to any changes in 
the PSF FHWM due to any differences in observing 
conditions. Furthermore, the PSF will only strongly 
affect the core region of the galaxy profile, which was 
not used in the final fit. 

A mask was needed to exclude other objects in 
the field affecting the profile fit. For the most 
of the field, objects were detected using SExtractor 
flBertin fc Arnouts 111996ft . However objects close to 
the core of NGC 4696 were not correctly detected, so 



objects in a region out to ~90 arcsec from the centre 
of the galaxy were detected and masked by eye. 

The galaxy surface brightness distribution was fit- 
ted using GALFIT. The models considered were a 
de Vaucouleurs profile and the more general Sersic pro- 
file. Different initial parameters were used to test the 
stability of the fitting algorithm. The 'sigma' image 
obtained from the combining process was used as the 
input error assigned to each pixel. 

There is a strong dependence of the fit results on 
the correct determination of the sky level. We there- 
fore kept the sky level as an additional free parameter 
in all the fits. In order to test the dependence on the 
sky level used, we also fitted profiles with the sky level 
held constant at ±0.2% of the best fitting value, to 
reflect the 0.4% variation seen across the image. The 
results are given in Table 1. 

The geometric parameters of the galaxy (axis ratio 
and orientation of the semi-major axis) were obtained 
from the initial fits. The axis ratio is b/a = 0.83±0.01 
and the position angle is 9 — 100° ±2° (measured north 
through east). These values were kept constant in the 
following fits, and used to construct masks limiting the 
region to fit, setting both the minimum and maximum 
radius for each fit. 

Masking the core of the galaxy was essential, as 
some complex structure was observed in that region 
that is unrelated to th e rest of the galaxy's profile 
(ICrawford et al. 1l2005T) . The region masked was an 
ellipse with axis ratio and orientation as above and 
semi-major axis of 18 arcsec. 

The goodness-of-fit can be evaluated from the x 2 
of the fit and the number of degrees of freedom (v) 
in the model. The value of \ 2 ls calculated directly 
by GALFIT. However, due to the correlation between 
neighbouring pixels caused by the PSF, v can not be 
calculated directly as the number of pixels minus the 
number of parameters. To overcome this, artificial im- 
ages were created from models with the same param- 
eters as the fits obtained and noise added at the same 
level as our data. These images were then smoothed 
using a PSF of the same FWHM as observed. Fits 
obtained from the resulting smoothed images give an 
estimate of the effect of the correlation. Values of the 
reduced x 2 (xt = X* l v ) quoted below take this cor- 
rection into account. 



5 Results and Discussion 

Successive fits were performed in regions limited by 
an ellipse with the same geometric parameters as the 
galaxy and a varying semi-major axis. The maximum 
distance at which the fit was reliable was estimated us- 
ing three indicators: the values of x2, visual inspection 
of the model-subtracted image, and the variation of 
flux with angle in successive elliptical rings. The max- 
imum allowed value for the semi-major axis distance 
was 180 arcsec (37 h^Q kpc) . All the results below refer 
only to fits within this region. 

The total magnitude (m), central surface bright- 
ness (fio), effective radius (r e ), and surface brightness 
at r e (fj, e ) obtained for de Vaucouleurs and Sersic pro- 
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Figure 2: Comparison between the observed pro- 
file and the fitted de Vaucouleurs model. Each 
point in the plot corresponds to the average sur- 
face brightness of the pixels within an elliptical 
annulus (obtained with IRAF task ellipse); r is 
the semi-major axis distance to that annulus. The 
error bars correspond to the dispersion in surface 
brightness of these pixels. The dashed (and dot- 
ted) lines correspond, respectively, to the values of 
r e and fi e and their uncertainties. 



files fitted to the studied region are shown in Tabled 
The extra free parameter in the Sersic model (the in- 
dex n) scarcely improves the value of \ 2 , and the index 
obtained is very close to the de Vaucouleurs value of 4. 
We therefore conclude that this extra free parameter 
is not needed, and that the light profile of NGC 4696 
is well described (in the region within 180 arcsec) by 
a de Vaucouleurs law. A plot of the observed surface 
brightness profile of NGC 4696, together with this best 
fit, is shown in Figure |5] 

The fit was repeated introducing variations of the 
geometric parameters as well as the masking of the 
core region in order to test the robustness of the re- 
sults. The main differences were found when varying 
the core masking, probably due to the presence of some 
extra structure close to the centre of the galaxy or un- 
certainties in the PSF. Errors quoted in TableQcorre- 
spond to these variations. The differences found when 
varying the geometric parameters (orientation and axis 
ratio) were smaller than these errors. 

The profile obtained was compared to that ob- 
tained in the if -band from 2MASS. The comparison 
is shown in Figure El There is a clear difference in the 
slopes of the profiles beyond the galaxy core. 2MASS 
fit a Sersic profile with an index of n = 1.11, whereas 
we find a profile that is well fit by an n — 4 de Vau- 
couleurs model. The difference between our result and 



Figure 3: Comparison with 2MASS surface bright- 
ness profile. The points and solid line are the 
2MASS profile and the dashed line is our fitted 
de Vaucouleurs profile (with parameters as in Ta- 
ble EJ). 



the 2MASS result is likely to be a consequence of the 
sky background estimation. It was found that the 
slope of the outer profile was very sensitive to the pre- 
cise sky level used. If we subtract a sky level higher by 
only « 0.3% than that obtained from the best fit, the 
difference was essentially eliminated. Note that this 
difference could be explained by the large scale varia- 
tions of « 0.4% observed in our final image, whereas 
2MASS has more optimally observed flat-fields. How- 
ever, our deeper observations are in better agreement 
with p revious work done in the optical; ISchombert I 
(1987) finds a good fi t with a de Vaucou l eurs model 
to r « 5 arcmin, and Uerien fc Dressier I (|l997f i fit a 
Sersic profile with n = 2.04. 

The extremely sensitive nature of the fits to the 
sky level (see Table 1) means that the previous un- 
certainty in the profile of NGC 4696 is still unresolved. 
Future observations could improve on the present work 
through either off-source blank field observations in- 
terleaved with the galaxy observations, or drift scan 
observations, as in 2MASS but deeper, to better esti- 
mate the flat-field and reduce the uncertainty in the 
sky levels. 



6 Conclusions 

The NIR light profile of NGC 4696 in the //-band was 
measured out to a distance of 180 arcsec (37 h^ kpc) 
along the semi-major axis. This limit was imposed 
by systematic errors in the flat-fielding process. The 
obtained profile is well described by a de Vaucouleurs 
law, with no signs of either truncation or an extended 
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Tabic 1: Surface Brightness Profile Parameters 



Parameter 


de Vaucoulcurs 


bersic 


bystematic errors 


m (mag) 


D.Zo =t U.UD 


^ i o _i_ n in 


(A Qfi O/l \ 
^4.oD, D.y4j 


/xo l-uicig aicsee j 


19 OQ^n 07 


11 q _i_ r\ c 
11. o zc u.o 




[i e (mag arcsec~ 2 ) 


20.62 ±0.07 


20.90 ±0.22 


(23.75, 19.36) 


n (Sersic index) 


4 


4.3 ±0.3 


(7.3, 2.4) 


r e (arcsec) 


172 ± 5 


204 ± 20 


(1199, 76) 


r e (h^ kpc) 


35.3 ± 1.0 


42 ±4 


(246, 16) 




1.170 


1.167 





"Limits based on fitting a Sersic profile with sky values fixed at ±0.2% the best fitting value. 



halo. These observations agree, up to an uncertainty of 
0.3% in the sky level, with the if-band profile obtained 
from 2MASS. 
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